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Abstract 

In this paper, we study the Cauchy problem of the linearized kinetic equations for the models of Marie and Anderson- 
Witting, and compare these dispersion relations with the 14-moment theory. First, we propose a modification of 
the Marie model to improve the resultant transport coefficients in accord with those obtained by the full Boltzmann 
equation. Using the modified Marie model and Anderson-Witting model, we calculate dispersion relations that are 
kinetically correct within the validity of the BGK approximation. The 14-moment theory that includes the time deriva- 
tive of dissipation currents has causal structure, in contrast to the acausal first-order Chapman-Enskog approximation. 
However, the dispersion relation of the 14-moment theory does not accurately describe the result of the kinetic equa- 
tion. Thus, our calculation indicates that keeping these second-order terms does not simply correspond to improving 
the physical description of the relativistic hydrodynamics. 
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1. INTRODUCTION 

Recently, the interest in relativistic dissipative flu- 
ids in astrophysics and nuclear physics has increased. 
Relativistic dissipative fluid equations have many fea- 
tures that do not appear in the case of nonrelativistic 
fluid. The most basic difference is the fact that in the 
presence of heat flux, the fluid velocity cannot be de- 
fined uniquely. There are two well-known definitions 
of fluid velocity; Eckart velocity UJ that is parallel to 
particle flow and Landau-Lifshitz velocity ||2[] that is 
parallel to energy flow. In addition, it is well known 
that standard first-order relativistic Navier-Stokes hy- 
drodynamics exhibits fatal problems regarding causality 
and stability, that is, small perturbations to the uniform 
static states grow exponentially JslQl- Currently, the 
most widely accepted and studied theory is the second- 
order Israel-Stewart (IS) approach f?) based on the 14- 
moment method |6]. Unfortunately, this theory is in- 
convenient for practical use because we have to re- 
store so many terms that are second-order in deviations 
from equilibrium, that is, the time derivative of dissi- 
pation terms and the products of gradients of dissipa- 



tive quantities. However, because of the recent finding 
of the strongly coupled quark-gluon plasma (sQGP) in 
the Relativistic Heavy-Ion Collider (RHIC), description 
by relativistic hydrodynamics equations have been vig- 
orously studied in the context of nuclear physics ||7|], 
and application of IS theory has just begun l'^. Re- 
cently, a new approach to relativistic dissipative fluid 
equation has been shown by Tsumura, Kunihiro, and 
Ohnishi nil . They use the renormalization- 

group method for obtaining fluid equation from the 
Boltzmann equation, and obtained equation is different 
from both Eckart and Landau-Lifshitz equation. 

Microscopic phenomena are accurately described by 
the Boltzmann equation. However, it is very diffi- 
cult to solve since its collision term depends on the 
product of the distribution functions. Consequently, 
a simpler approximation for the collision term has 
been proposed; the most widely used relativistic ki- 
netic model equations are those of Marie |.12.] and 
Anderson- Witting ifisll . The Marie model is an ex- 
tension of the nonrelativistic Bhatnagar-Gross-Krook 
(BGK) model 1 14] to the relativistic case and is de- 
scribed in the Eckart frame 1 1 ] . The Anderson- Witting 
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a boundary value problem by Cercignani and Majo- 
rana 241. To understand the dynamics as a Cauchy 



model is another extension and is described in the 
Landau-Lifshitz frame |2]. Of the two, the Anderson- 
Witting model is widely used ifisl Eoll because the 
Marie model has undesirable properties; for example, 
the transport coefficients obtained by the Marie model 
do not agree with those obtained by the full Boltzmann 
equation [6]. 

In this paper, we compare the dynamics described by 
the 14-moment theory with that of the kinetic model 
equation and test the applicability of the IS approach. 
To do numerical simulation of relativistic dissipative 
fluid, we should know how to treat the small second- 
order terms and how to determine appropriate values of 
new coefficients, which urges us to check how impor- 
tant these terms are. To make the problem tractable, 
we study linear perturbation and compare the solutions 
of the dispersion relation. The dispersion relations of 
the relativistic kinetic equations have been studied as 
mndary \ 

problem, we solve the dispersion relations with respect 
to cl>. In addition, we modify the problematic properties 
of the Marie model and use the modified model equation 
to analyze the Eckart description. 

This paper is organized as follows: in Sec. |2] we 
introduce the kinetic models of Marie and Anderson- 
Witting. Then, we modify the Marie model and obtain 
the dispersion relations. In Sec. [3] we solve the disper- 
sion relations numerically with respect to o) and present 
our results. In Sec. ID we discuss the properties of the 
Marie and Anderson- Witting models. In addition, we 
analyze the asymptotic behavior of the dispersion rela- 
tions. First, we study the long wavelength limit and then 
we study the short wavelength and high frequency lim- 
its. We solve the dispersion relation of the 14-moment 
theory and compare it with the dispersion relations of 
the kinetic model equation. 

2. THE LINEARIZED KINETIC EQUATIONS 
AND THE DISPERSION RELATIONS 

In this section, we derive the dispersion relations of 
the relativistic kinetic models of Marie and Anderson- 
Witting. Throughout this paper, we use the units 



(1) 



where c is the velocity of light, and ks the Boltzmann 
constant. 

In Cartesian coordinates, the Minkowski metric ten- 
sor rifjy is given by 



Variables indicated by Greek letters take values from 
to 3, and those indicated by Roman letters take values 
from 1 to 3. 

2.1. MODIFICATION OF THE BGK MODEL OF 
MARLE 



Marie Ill2ll has proposed the following form of the 
kinetic model equation. 



(/(f,x,p)-^(f,x,p)) 



(3) 



where tm is a characteristic time on the order of the 
mean flight time (see below for its physical interpreta- 
tion), m is the rest mass of a particle of the relativistic 
gas, and feq is the local equilibrium distribution func- 
tion. 

Using Eq. Q, we obtain the following form of the 
kinetic equation 



{f-.fe,\ (4) 

(5) 



The Marie model is an extension of the nonrelativis- 
tic BGK model to the relativistic case. The transport 
coefficients for the Marie model equation reproduce the 
nonrelativistic results in the limiting case of low tem- 
perature. It is, however, well known that in the limit- 
ing case of high temperature, the transport coefficients 
of the Marie model differ from those found for hard- 
sphere particles obtained by the full Boltzmann equa- 
tion |6]. More precisely, if we express the transport co- 
efficients (oc Tm) as a function of ^ = mjT, the transport 
coefficients of the Marie model behave as 1 of those 
found by the Boltzmann equation for hard-sphere parti- 
cles in the limit of high temperature. For this problem, 
we should recall that the transport coefficients are gen- 
erally proportional to the relaxation time tm, and Eq. Q 
contains tm as a parameter of the BGK model. This in- 
dicates that the appropriate value of tm is different from 
the physical relaxation timescale Treiax by a factor that 
becomes unity in the low temperature limit and becomes 
^ in the high temperature limit. We discuss this new in- 
terpretation of Tm. 

First, we clarify the meaning of the parameter t in 
the BGK model. In the nonrelativistic BGK model, the 
parameter r is equivalent to the relaxation time. The 
nonrelativistic kinetic equation of the BGK model is 



77^v = diag(l, -1,-1,-1). 



(2) 



B \ 1 

- + V ■ V / = --(/(f, X, V) - X, V)). 



(6) 
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If the one-particle distribution function / does not de- 
pend on the spatial coordinates, Eq. (|6|l reduces to the 
ordinary first-order differential equation, and we can ob- 
tain the formal solution 



fit) 



m + - f e''l'fe,(t')dt' 

T Jo 



(7) 



This equation indicates that t is the relaxation time of 
the distribution function. 

Next, we consider the relativistic BGK model of 
Marie. The kinetic equation of the Marie model is 



--^(/(f,x,v)-/,,(f,x,v)). (8) 

If we assume that the one-particle distribution func- 
tion / does not depend on the spatial coordinates, the 
formal solution of Eq. ([8]l is 



/(f) = 



/(0) + 



-Tm- 



j_ r //TM, 

Tm* Jo 



feq(t')dt' 



-t/Tu, 



(9) 



(10) 



This indicates that in a general inertial frame, the relax- 
ation time is not tm but tm*, and tm is the relaxation 
time in the rest frame where the momentum of particles 
is p = 0. More precisely, if we employ the particle's rest 
frame where p = 0, Eq. (|8]l becomes 



d m 

-fit, X, 0) = (fit, X, 0) - fe.it, X, 0)). 

dt Tm 



(11) 



This is the same equation as in the nonrelativistic BGK 
model, indicating that only in this frame does tm be- 
come the relaxation time. 

Although the transport coefficients of the Marie 
model are expressed in a form proportional to tm in 
the literature Jid, the above explanation shows that 
we should use tm* as the relaxation time instead of 
Tm- However, tm* depends on the momentum so 
Tm* cannot appear in macroscopic descriptions, such as 
transport coefficients. For this reason, we have to con- 
sider the true relaxation time Treia.x, to which the trans- 
port coefficients should be proportional, and relate it 
to the BGK parameter of the Marie model tm- The 
above discussion suggests that we may regard 1/Treiax 
as {1/tm*}, and we can consider Tieiax as the effective 
relaxation time in general frames. Using the local equi- 
librium distribution function, tm is 



Tm 



~ I n Jeq'^relax ~ , ^\ "^relaxi 



(12) 



where K„ is the second kind modified Bessel function 
of order n. The correction K\{^) I K2{0 becomes 1 in 
the limit of large ( and (jl when f is nearly 0. This 
indicates that this function has the desired properties. In 
the following, we use this tm as the BGK parameter of 
the Marie model. 

In above discussion, we assume that the physical sys- 
tem is not far from equilibrium state, and calculate the 
average of tm with respect to the local equilibrium dis- 
tribution function f^q. Though this cannot give the cor- 
rect Tm in the general case, it is a good approximation 
for linear perturbation about the local equilibrium dis- 
tribution function. 

2.2. THE LINEARIZED KINETIC EQUATION AND 
DISPERSION RELATION OF THE MODIFIED 
MARLE MODEL 
In this section, we derive the dispersion relation of the 
modified relativistic kinetic model of Marie. To obtain 
the dispersion relation, we apply an approach similar to 
that in the work of Cercignani and Majorana [21]. 

When there is no external field, the equation of the 
modified Marie model is given by 



D m 
Ds Tm 



(13) 



= = — V)/. (14) 

In Eq. ( fTlT l. Tm is the relaxation time modified in 
Sec. 12.11 and f^g represents the local Maxwell- Jiittner 
distribution function 



/c^(f,x,p) = 



n(t, x) 



X exp 



47Tm^T(t,x)K2(^(t,x)) 

PfiUf'it, x) 



T(t,x) 



(15) 



(16) 



where m is the mass of the particle, and T is the temper- 
ature. 

Eq. (HJI ) is a nonlinear equation for f{t, x, p) because 
of the nonlinear dependence of f^g on / through the fol- 
lowing conditions called the matching conditions: 



fe,-f)/-^-0. 



(17) 
(18) 



To obtain the dispersion relation, we start by expand- 
ing the distribution function around a global equilibrium 
state /o(p), 



Sf ^ f - fo, 6feq = feq - fo- 



(19) 
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Then, the linearized kinetic equation of the modified 
Marie model is given by 



tm* 



tm* - — Tm- 
m 

We assume a solution in the following form: 



6f = dfe-""'' = 5fe 
Then, Eq. (l20b reduces to 



— icj + ik ■ V ) 5/ = —dfeq. 

Tm* / Tm* 



(20) 
(21) 

(22) 
(23) 



We consider an equilibrium background state in which 
the fluid is at rest, so that m'' - (1,0) and du'' - (0, 6u) 
owing to the relation m'^^m^ = 0. Then, Sfeq is given by 



6n 
n 



-1 



+ 

T K2 ] T 



p ■ dn 



/() 



■ exp 



' T 



(24) 
(25) 



where Sn is the space component of the Eckart velocity. 



as explained in Sec. Appendix B and K'^^ is the deriva- 
tive of Kn with respect to f . 

Using the matching conditions, we can rewrite 6p, 6vl, 
and ST as the integrals of 6f: 



5n(f,x) = J ^A/, 

5u(f,x) = — -^^Sf, 
n J 



6T{t,x) 

J 1 
Eq. (l23T l becomes 



(26) 
(27) 
(28) 



/ w + /k ■ V 5 /(p) 



(29) 



^^V^' /o(P) 



P P 



«(i-/:i^//:2)(3 + ^2 + ^j^/^2) 

\ T K2^ \ T K2^ 



In the following, we take T,-eiax as a unit of time: 

CJTrelax ^ W, T^e/aj:^ ^ ^. (30) 

Finally, the linearized equation of the BGK model of 
Marie is 



/o(P) 

l-(/c.-/k.i)|| 



K{p,p') = 



P P' 



(31) 
(32) 



T 

+ — 



n(l-KiaK2)(3 + e + KiaK2) 



where z = p**/?^. This equation make sense only when 
1 - (iLo - ik ■ ^) ^ 7^ 0; we explain the case where 

1 - [ioj - /k ■ §f = later 

Eq. dJTT l is the homogeneous Fredholm integral equa- 
tion of the second kind. In particular, the kernel function 
K{p, p') can be separated with respect to the variables p 
and p'. Thus, this equation can be solved according to a 
general procedure. 

First, we integrate Eq. ( |3T] ) with respect to p. Then, 
we multiply by Kik/( and the equation reduces to 



6n ST 
/ii— +/i2k-5u + /i3— =0, 
n T 



where 



/ii = K^k - (K2P(Q) - nK2e''^, 



(33) 

(34) 
(35) 



■ nK2e-''j^^-iu>(\+d) + ^^C 



/i3 = (-3 - ^^j {kK\ - (K2P(Q)) (36) 



+ ^K2 (kK, - (P{ 1)) - 7TK2e-d \d-2-j-n 



P(n) is 



Pin) 
b 



dy e ^-'y" arctan 



and d is 

d^ — 

d^(, 



k\K, 

K2( 

/:ilm(w) 



luy 



if Im(w) <- — , 

if Im(w) >~. 

^1 



(37) 
(38) 

(39) 
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The derivation of the correction term c is explained in 



Sec. Appendix A. 2 



Next, we muhiply Eq. ( |3TI ) by k ■ p = kp'^ and inte- 
grate with respect to p. Then, we multiply by ^Ki and 
the equation reduces to 

6n ST 
/21— +/22k-(5u + /23— = 0, (40) 
n T 



where 



(41) 



-d 



w(l +d) + i—i}' 



/.4{^.(l-M-g(|p(0)-/.P(l)) (42) 



+ -e-''^\(J(cf + 2d + 2) 
k k^ ^ 



(43) 



K\k \ K2 
+ -^(-3P(l)-^P(l) + fP(2) 

{K' 
(l+d)—( + d^+d+l 
K2 

K2J, K' 

^'kA'^T2^ 

Next, we multiply Eq. dSTT l by k x p s kp^ and inte- 
grate with respect to p. Then we multiply by K\k and 
the equation reduces to 



h±Su^ - 0, 



where 

/i_L = 2kKi 



(44) 



(45) 



{^j^^V- l)}arctan 



2iK2 
k^Ki 



-^^e-'{id'+2d + 2)ll-^ 



(Ci)(l + d) + 



Eq. (l40l i corresponds to the longitudinal mode (6ujc + 
0, 5m J. = 0), and Eq. ( l44l l corresponds to the transverse 
mode (5mv = 0, 6u^ + 0). 

Finally, we multiply Eq. dSTT i by /^"^ and integrate with 
respect to p. We multiply by K\k and the equation re- 
duces to 



where 
/31 

/32 



6n 5T 

h\— + hi^- + hi- 
ll T 



■ (P{l)-K,k+-e-\\+d), 



0, 



P(l)-iuP(2) 
K2 



+ -e-"- \ -icjid' +2d + 2) + —^(1 + d) 



(46) 

(47) 
(48) 

(49) 



/33 =^ 1^-3^(1) -■^mi) + m2) 

+ -e-'' ((1 + d)—C + d^ +d+l 
(, \ K2 



If the determinant of the above homogeneous system 
is set to equal zero, the following dispersion relation is 



obtained: 






^11 hi hi 







hi hi hi 







hi hi hi 


















/x± 


This condition implies either 




h\ hi hi 






hi hi hi 


= 0, 




hi hi hi 





(50) 



= 0. 



(51) 



(52) 



Using this dispersion relation, we can obtain 6f in the 
form 



C„/o(v) 



(53) 



n \ T K2^ T 



P 6Ua 



,-;(w„/+k-x) 



where C„ is a constant coefficient, and 6n^^^,6Tij^^, and 
dUai,, are eigenfunctions obtained from the dispersion re- 
lations. 
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If 1 Jtm* — + /k ■ V = 0, the mode becomes contin- 
uous II22II23II . According to Eq. ( l23T l. the eigenfunction 
for this mode satisfies the equation 



Then, Eq. ( 1581 ) reduces to 



= 6f,, 



(54) 



This mode represents the decay of the moments of / 
with vanishing 6feq, i.e., 6n - ST - 0, 5u = 0. 

Unhke the case of the nonrelativistic BGK model and 
Anderson-Witting model, the decay rate of this contin- 
uous spectrum is not constant but depends on 112 ill . 



i -/w + /k-vj5/= i^^. (60) 



We consider an equilibrium background state, in 
which the fluid is at rest, so that u'^ - (1,0) and du^^ - 
(0, Su) due to the relation m''i5m^ = 0. Then, 6feq is 



n \ T K2 j T T 



■ exp 



2.3. THE LINEARIZED KINETIC EQUATION AND 
DISPERSION RELATION OF THE ANDERSON- 
WITTING MODEL 

In this section, we derive the dispersion relation of the 
relativistic kinetic model of Anderson- Witting. To ob- 
tain the dispersion relation, we apply an approach simi- 
lar to that in the work of Cercignani and Majorana 12411 . 

When there is no external field, the equation of the 
Anderson- Witting model [ 1 3i1 is given by 



, (61) 



(62) 



Ds ^ 



(55) 



-UyP 



vf fee 



47Tm^TK2(0 

where 6u is the space component of the Landau-Lifshitz 
velocity, as explained in Sec. [Appendix B| 

Using the matching conditions, we can rewrite 6p, 6u, 
and ST as integrals of Sf: 

6n(t,x)^ J ^p"6f, (63) 

6n(t,x)^-- [ ^pSf, (64) 

-^P° (65) 
-l+K'2aK2+pVT 



where feg is the local equilibrium distribution function 
defined by Eq. ( fTSl l. 

As in the Marie model, Eq. {55[ is a nonlinear equa- 
tion for /(f , X, p) because of the nonlinear dependence 
of ft,q on / through the following matching conditions: 



Then, Eq. ( fSSl ) becomes 



iaj + ik ■\]6f(p) 



(66) 



-I 



d'p' /o(p) 



P P 



n^3 +i P 

«v -^py [feq - /) = 0, (56) n{\-K,aK2){-i+e + Ki (IK2) 

I rfi K' \ I K' 



where ij/ is defined by Eq. ( fTST l. 

To obtain the dispersion relation, we start by expand- 
ing the distribution function around a global equilibrium 
state /o(p). 



T K 



T K2 

In the following, we take r as a unit of time: 
(jjT Lo, rk ^ k. 



(67) 



Sf ^ f - fo, Sfeq = feq ' /() 



(57) 



Finally, the linearized equation of the BGK model of 
Anderson- Witting is 



The kinetic equation of the Anderson- Witting model 
Eq. ( fSST l reduces to 



5/(p) = J ^/:(p,p')W), 



d \ Sf- 6 fee 



in a linear approximation. 

We assume a solution of the form 



(58) 



^(P,P') = 



/o(p) 



I - io) + ik ■ 



P P 



n 

„/0 



(68) 
(69) 



n(l-KiaK2)i3+^+KiaK2) 



Sf = ^/e-'*"-'"' = Sf 



(59) 



Sfip'l 
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where z - IT . This equation make sense only when 
1 - + /k ■ 4 ^ . When 1 - /w + /k ■ 4 = , the 
mode becomes continuous, as in the Marie model. 

Eq. (|68] | is a homogeneous Fredholm integral equa- 
tion of the second kind. In particular, the kernel function 
K{p, p') can be separated with respect to the variables p 
and p', and this equation can be solved according to a 
general procedure. 

First, we multiply Eq. ( |68] | by and integrate with 
respect to p. Then, we multiply by and the equation 
reduces to 



5n 6T 
/ii— +/i2k-(5u + /i3— =0, 
n T 

where 



and Qin) and b are defined as follows: 



b = 



dy e ^^'y" arctan 
1 - ia> 



by 



(74) 
(75) 



Next, we multiply Eq. ( |68] | by k ■ p and integrate 
with respect to p. Then, we multiply by K2, and the 
equation reduces to 



6n 6T 

hi — +/22k-5u + /23 — 

n T 



0, 



where 



/21 = 3K2 + (K, - bfQ{3), 
b 

hi = 
+ 3^K, 



(Ki - -[{12 + C'')K2 

^r'e(4)}], 



/23 =-^^1^(3 + ^1 + (3 + ^2)^2 
■be\-3Q0)-j-^^Q0) + m^) 



(76) 

(77) 
(78) 

(79) 



Next, we multiply Eq. ( |68] ) by k x p = kp^p-^ and 
integrate with respect to p. Then, we multiply by 2kK2, 
and the equation reduces to 



where 

/^^ = 2k^K, - [(b^ + 1)2(4) - 2(2)} (81) 
+ fe {(12 + ^2)^:2 + 3^^:2}, 

Finally, we multiply Eq. ( |68] | by (/?*')2 and integrate 
with respect to p. Then we multiply by K2k, and the 
equation reduces to 



(70) 


where 




/31 


(71) 


/32 


(72) 


^33 


(73) 





Sn ST 
hi—+ /32k ■ 6u + /33— = 0, 
n T 



k{3K2+(Ki)-fQ{3), 
-^[b^'Q(4)- (^^ + l2)K2-3iKi]., 



K2 



(82) 

(83) 
(84) 

(85) 



■r 1-32(3) --^2(3)^ + ^2(4) 



If the determinant of the above homogeneous system 
is set to zero, we can obtain dispersion relation the same 
as Eqs. (|50]l, (|5B, and (|52]i. 

Using this dispersion relation, we can obtain Sf in the 
form 



sm = 2 — 



C„/o(v) 



ia> + ;k ■ -ir 



(86) 



n \ T K2 



6T„. 



hi.5u±_ - 0, 



(80) 



where C„ is a constant coefficient, and 6n^^^,6Ti^^^, and 
5Uw,_ are eigenfunctions obtained from the dispersion re- 
lations. 

3. RESULTS 

3.1. MARLE MODEL 

In this section, we show the dispersion relations of 
the modified Marie model obtained in the previous sec- 
tions. We solve the dispersion relations numerically; the 
results are shown below. First, we show the thermal 
conduction mode in Figs.[Tl|2] and [3] At long wave- 
lengths, the decay rate is proportional to k^ and repro- 
duces the result obtained by the first-order Chapman- 
Enskog approximation. Note that the decay rate of the 
thermal conduction mode diverges at finite wavelengths 
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0.01 0.1 1 10 100 

k T c 



Figure 1: Decay rate of the thermal conduction mode in the non- 
relativistic case: mIT = 100. The black curve corresponds to the 
Anderson-Witting model while the thick gray curve corresponds to 
the modified Marie model. 




0.01 0.1 1 10 

k T c 



Figure 2: Decay rate of the thermal conduction mode of the modified 
Marie model in the relativistic case: mjT = 1. The gray zone repre- 
sents the region in which the BGK approximation is expected to be 
worse; the black zone represents the region in which we cannot use 
the BGK approximation. 

in the relativistic cases. This may be equivalent to the 
critical frequency of the thermal wave mode predicted 
in the work of Cercignani and Majorana [24]. We will 
return to this problem later. Second, we show the sound 
wave mode in Figs. 21 15] and |6] As in the thermal 
conduction mode, at the long wavelengths the decay 
rate is proportional to k^. In the relativistic and ultra- 
relativistic cases, the phase velocity becomes larger than 
light velocity at some wavelength, and we stop the cal- 
culation because the physical collision term produces a 
phase speed less than light velocity [28.1 . As in the case 
of thermal conduction mode, this may be equivalent to 
the critical frequenc y pr edicted in the work of Cercig- 
nani and Majorana 12411 . Numerically we obtain that 
kmax - 80 when ^ = 5 and k,„ax > 100 when ^ = 10 



100 




k X c 

Figure 3: Decay rate of the thermal conduction mode of the modified 
Marie model in the ultra-relativistic case: m/T = 0.01. 
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0.01 0.1 1 10 100 

k T c 

Figure 4: Dispersion relation of the sound wave mode in the non- 
relativistic case: mjT = 100. Solid curve represents the decay rate 
(-Im tS); dotted curve represents the frequency Re u). The black curve 
corresponds to the Anderson- Witting model while the thick gray curve 
corresponds to the modified Marie model. 



where k is the maximum wavelength of applicability of 
BGK model. 

Finally, we show the shear flow mode in Figs. [T] [8] 
and|9] As in the thermal conduction mode, at the long 
wavelengths the decay rate is proportional to k^. The 
dispersion relation for shear flow has only a decay rate, 
indicating that in rarefied gas shear flow cannot propa- 
gate. 

In nonrelativistic case, the relevancy of the adopted 
equation can be checked by comparing its dispersion 
relations to experimental data of attenuation rate and 
phase velocity of sound wave. Unfortunately, corre- 
sponding experiments are very difficult in relativistic 
regime, and we cannot compare our results to exper- 
imental data. However, our results of nonrelativistic 
case - 100) reproduce the dispersion relation of non- 
relativistic BGK lll4l E3I1 that agrees with experimen- 
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Figure 5: Dispersion relation of the sound wave mode of the modified 
Marie model in the relativistic case: m/T = 1. Solid curve represents 
the decay rate (-Im cj); dotted curve represents the frequency Re lo. 
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Figure 6: Dispersion relation of the sound wave mode of the mod- 
ified Marie model in the ultra-relativistic case: m/T = 0.01. Solid 
curve represents the decay rate (-Im oj); dotted curve represents the 
frequency Re at. 

tal data |27] even in short- wavelength regime. For this 
reason, we expect that our relativistic dispersion rela- 
tions should be correct even in relativistic regime, at 
least qualitatively. 

3.2. ANDERSON AND WITTING MODEL 

In this section, we show the dispersion relations of 
the Anderson- Witting model obtained in previous sec- 
tions. First, we show the thermal conduction mode in 
Figs. [To] and [TT] The nonrelativistic limit m/T - 100 is 
given in Fig. [1] At long wavelengths, the decay rate is 
proportional to and reproduces the result obtained by 
the first-order Chapman-Enskog approximation. As in 
the Marie model, the decay rate diverges at finite wave- 
lengths. 

Second, we show the sound wave mode in Figs. [12] 
and [13] The nonrelativistic limit m/T = 100 is given in 
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Figure 7: Decay rate of the shear flow mode in the nonrelativistic 
case: m/T = 100. The black curve coiresponds to the Anderson- 
Witting model while the thick gray curve corresponds to the modified 
Marie model. 
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Figure 8: Decay rate of the shear flow mode of the modified Marie 
model in the relativistic case: in/T = 1. 

Fig. [4] As in the thermal conduction mode, at the long 
wavelengths the decay rate is proportional to k^. In the 
relativistic and ultra-relativistic case, the phase velocity 
becomes larger than light velocity at some wavelength, 
and we stop the calculation as in the Marie's model. In 
this case, the wave number k^ax at which phase velocity 
becomes faster than light is ki„„x ^ 64 when ^ = 5 and 
kmax > 100 when f = 10. 

Finally, we show the shear flow mode in Figs. [T4l and 
[TS] The nonrelativistic limit m/T = 100 is given in Fig. 
[T] As in the thermal conduction mode, at the long 
wavelengths the decay rate is proportional to k^. The 
dispersion relation for shear flow has only a decay rate, 
indicating that in rarefied gas shear flow cannot propa- 
gate. 

In the Anderson-Witting model, we find the kinetic 
mode. Figs. [16] and [T7] show the transverse kinetic 
mode, and Figs.[T8]and[T9]show the longitudinal kinetic 



9 




le-05 I — — — — I 

0.01 0.1 1 10 100 

k T c 

Figure 9: Decay rate of the shear flow mode of the modified Marie 
model in the ultra-relativistic case: mjT = 0.01. 
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Figure 10: Decay rate of the thermal conduction mode of the 
Anderson-Witting model in the relativistic case: mjT = 1. 

mode. Due to the numerical difficulties, we can- 
not find the kinetic mode in the ultra-relativistic case: 
because the phase velocity of the longitudinal kinetic 
mode in the relativistic case is a little faster than that 
of light, this mode is physically incorrect. We expect 
that this problem results from the accuracy of the BGK 
model and does not show that there is any longitudinal 
kinetic mode of in the relativistic case. 

As for the Marie model, we cannot compare these 
dispersion relations and experiment; however, we ex- 
pect that our relativistic dispersion relations is correct 
even in relativistic regime, at least qualitatively for the 
same reason of Marie model. 

4. DISCUSSION 

4.1. ANALYSIS IN THE SHORT WAVE LENGTH AND 
HIGH FREQUENCY LIMIT 

In this section, we analyze the dispersion relation 
for the short wavelength limit of the shear flow mode 
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Figure 11: Decay rate of the thermal conduction mode of the 
Anderson- Witting model in the ultra-relativistic case: m/T = 0.01. 
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Figure 12: Dispersion relation of the sound wave mode of the 
Anderson- Witting model in the relativistic case: m/T = 1. Solid 
curve represents the decay rate (-Im oj); dotted curve represents the 
frequency Re w. 



and the high frequency limit of the thermal conduction 
mode. In particular, divergence of the decay rate of the 
thermal conduction mode is not observed in the analysis 
of the nonrelativistic BGK model, so we have to check 
this divergence analytically. For simplicity we analyze 
the Anderson- Witting model. 



First, we analyze the dispersion relation for large 
wave number in the shear flow mode. From Sec. 12.31 
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Figure 13: Dispersion relation of the sound wave mode of the 
Anderson-Witting model in the ultra-relativistic case: m/T = 0.01. 
Solid curve represents the decay rate (-Im oj); dotted curve represents 
the frequency Re at. 
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Figure 14: Decay rate of the shear flow mode of the Anderson- Witting 
model in the relativistic case: m/T = 1. 
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Figure 15: Decay rate of the shear flow mode of the Anderson- Witting 
model in the ultra-relativistic case: m/T = 0.01. 
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Figure 16: Dispersion relation of the transverse kinetic mode of the 
Anderson- Witting model in the nonrelativistic case: m/T = 100. 
Solid curve represents the decay rate (-Im at); dotted curve represents 
the frequency Re at. 



the dispersion relation of the shear flow mode is 

n6u^^ J ^P±Sf (87) 

j p" 



feq 

\ — io) + ikVx 



\5n I p° K^\6T 
x\ — + 1-1 + — + —n — 
[n \ T K2 T 



(/y/l+h'- 



2n r 

I 



d^P P±feq 

1 - iu + ikv; 



■p±6uj 



b = 



'a 

1 — io) 



(88) 



In the above equation, we add correction terms of an- 



alytical continuation, because the decay rate -Im to is 
larger than 1/t in the short wavelength limit. 

In the high frequency limit, we use the following ap- 
proximations: 

I - lai + ikvx i<jL>\ 10) 



10) 



(90) 
(91) 



Using the above approximations and neglecting the 
terms higher than the third order of we reduce 

Eq. dHlll to 



1 : TT^a - 0, 



k fc2 



dz Z^fec, 



(92) 
(93) 
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Figure 17: Dispersion relation of tlie transverse kinetic mode of the 
Anderson- Witting model in the relativistic case: m/T = 1. Solid 
curve represents the decay rate (-Im at); dotted curve represents the 
frequency Re at. 



Figure 19: Dispersion relation of the longitudinal kinetic mode of 
the Anderson-Witting model in the relativistic case: m/T = 1. Solid 
curve represents the decay rate (-Im at); dotted curve represents the 
frequency Re oj. 
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Figure 18: Dispersion relation of the longitudinal kinetic mode of 
the Anderson- Witting model in the nonrelativistic case: m/T = 100. 
Solid curve represents the decay rate (-Im cj); dotted curve represents 
the frequency Re a). 



ber Eq. ^ is 



6n r t/ p 



feq 



iojikvx 



(95) 



T 5n 1 T 

{ct +2d + 2)— + ib{cF + 3ct + 6d)Su^ 
n 



ST 

+ + 3cP- + 6d + 6\ — 



d = 



(96) 



Finally, we obtain the dispersion relation of the shear 
flow mode in the high frequency limit, 



0) - -i ^Inak^^'^, 



(94) 



where we take the sign representing the decaying mode. 
This reproduces the results in Sec. 13.21 In the nonrela- 
tivistic regime, we cannot neglect the w/fc dependence 
in a, and the behavior of the shear flow mode in the 
large wave number limit becomes different from k^l^. 

Next, we consider the high frequency limit of the 
thermal conduction mode. Since the decay rate di- 
verges, we study only the highest-order terms in d. 



From Sec. 12.31 the conservation law for particle num- 



In the above equation, we add the correction terms of 
analytical continuation, because the decay rate -Im to 
is larger than 1/t in the large wave number limit. 

As in the shear flow mode, we expand the integrand 
in powers of u and neglect terms higher than second 
order on the right-hand side. Then, the above equation 
reduces to 



-1 + 



In \ 5n 

6mb 



(97) 



Kik^- 



Un K, ST 



K2keK2^'T 



where we consider the relativistic limit f <s; 1, so that 
we approximate exp[-(i] ^ 1. 

Similarly, the conservation of energy Eq. (l82b reduces 
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(98) 
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From Sec. Appendix B the dispersion relation of the 
Anderson-Witting model includes the conservation of 
energy Eq. (IB. 9b . so we use the conservation of energy 
instead of the conservation of momentum Eq. ( l76l l. 



0.01 
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0.0001 
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• iLo(c^,T6n + nCydT) + inhkdux — 0, 



(99) 



From the above equations, we can obtain the dispersion 
relation in the form 



.^^[lAcykKlne (100) 

- Ucyii^ {AK2 + IKi^i-l + 80)] + 0(a)) = 0, 

In the high frequency limit, terms of lower order than 
up' can be neglected. For the left-hand side to vanish, 
the coefficient of uP should be 0. Then we obtain 



'\K2 - ^Ki^ + lAKie 

iKie 



(101) 



If we insert ( - 0.01, we obtain k ^ 3.141 ■ ■ ■ . This re- 
produces the critical wave number for thermal conduc- 
tion in the Anderson-Witting model accurately. Simi- 
larly, if we insert ( - 100, we obtain k ~ lO'*^. This in- 
dicates that in the nonrelativistic regime the Anderson- 
Witting model does not effectively yield the critical 
wave number In deriving the above equation, we as- 
sume exp[-c/] ~ 1, so we can not reproduce the critical 
wave number of ^ - \ very well. More accurate anal- 
ysis reproduces the critical wave number of ^ = 1 as 
k ^ 4.958. 



4.2. COMPARISON TO 14-MOMENT EXPANSION 

In this section, we analyze the dispersion relation of 
the 14-moment theory and compare it with that of the 
BGK model. We assume that the relativistic gas is at 
rest, and we consider only longitudinal waves. In this 
case, the dispersion relation is given in the work of Cer- 
cignani and Kremer|'6','l6']. The IS equation is based on 
14-moment theory; thus, results obtained in this section 
can be applied to IS equation as well. 



Figure 20: Decay rate of the thermal conduction and kinetic modes of 
the 14-moment theory in the nonrelativistic case: m/T = 100. 



In the nonrelativistic limiting case ^ » 1 , the disper- 
sion relation reduces to 
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(107) 



We are interested in the Cauchy problem, so we solve 
the above equation with respect to co. The results are 
illustrated in Figs.|20]and|2T|in the case of ^ = m/T - 
100. Fig. |20] shows the heat conduction mode and its 
accompanying kinetic mode. This figure shows that at 
the long wavelengths, the decay rate of the heat conduc- 
tion mode of the 14-moment theory is proportional to 
k^ and reproduces the result of the first-order Chapman- 
Enskog approximation. At short wavelengths, the de- 
cay rate of the heat conduction mode has an upper limit 
and approaches the limit asymptotically. In addition, 
14-moment theory reproduces the kinetic mode. Fig. 1211 
shows that the decay rate of the sound wave mode has 
similar features. In comparison to the kinetic model 
equation and experiment |27], we find that in the nonrel- 
ativistic limit, the behavior of the decay rate of the fluid 
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Figure 21: Dispersion relation of the sound wave mode of the 14- 
moment theory in the nonrelativistic case: mjT = 100. Solid curve 
represents the decay rate (-Im tj); dotted curve represents the fre- 
quency Re 0). 



mode of 14-moment theory is consistent with the kinetic 
equation at long wavelengths but inconsistent at short 
wavelengths. In contrast to 14-moment theory, the BGK 
approximation reproduces the result of experiment 
qualitatively. In addition, the kinetic mode obtained by 
14-moment theory decrease with k in contrast to kinetic 
modes of BGK equation. 

Next, we consider the ultra-relativistic limit ^ <s; 1 . 
In this case, the dispersion relation reduces to 



225 , 35 2 5 

lU Cl> lO) 

16 4 4 



(108) 



( kc, 
[ u 
125, 



175, 



145 



25. 1 

— lOJ + 1 



.... 3 145 2 73. , ^ 

+ l(x) U lUJ -1-1=0. 

16 12 12 

The results are illustrated in Fig. |22] in the case of 
f - mjT - 0.01. Fig. |22] shows the heat conduc- 
tion mode and its accompanying kinetic mode. We 
could not calculate the sound wave mode very accu- 
rately because of the numerical difficulty in the complex 
Newton-Raphson method, but the behavior of the decay 
rate of the sound wave mode seems to be similar to that 
in the nonrelativistic case. These results indicate that 
at short wavelengths, the dispersion relation of the 14- 
moment theory is qualitatively diff^erent from that of the 
kinetic equation in the ultra-relativistic limit. 

In conclusion, the 14-moment theory is better than 
the first-order Chapman-Enskog approximation in the 
sense that this theory is causal and can describe the 
kinetic mode. However, the 14-moment theory can- 
not reproduce the result of the kinetic equation at short 
wavelengths or high frequencies, even in the nonrela- 
tivistic limit in contrast to the kinetic model equations. 



Figure 22: Decay rate of the thermal conduction and kinetic modes of 
the 14-moment theory in the ultra-relativistic case: m/T = 0.01. 



This indicates that the second-order dissipation terms 
do not reproduce kinetically correct results, and keep- 
ing these second-order terms may not necessarily im- 
prove the physical description of the fluid phenomena; 
but just make the mathematical form of fluid equations 
hyperbolic. 



5. CONCLUSION 

In this paper, we have solved the dispersion relation 
of the kinetic equations of the Marie and Anderson- 
Witting models with respect to w as a function of k 
since we are interested in the Cauchy problem. To ob- 
tain the dispersion relation, an approach similar to that 
in the work of Cercignani and Majorana f2T,'27l is ap- 
plied. To obtain an acceptable dispersion relation, we 
have modified the Marie model since it cannot repro- 
duce correct results in the relativistic case. Our disper- 
sion relation indicates that both kinetic model equations 
have a critical wavelength for the sound wave and ther- 
mal conduction modes; for the sound wave mode, the 
phase velocity exceeds the speed of light at that wave- 
length 1I28I1 . and for the thermal conduction mode, the 
decay rate diverges at that wavelength. 

We have solved the dispersion relation of the 14- 
moment theory with respect to w as a function 
of k and compared it with that of the kinetic model 
equations. The results show that the 14-moment theory 
reproduces the first-order Chapman-Enskog approxima- 
tion in the long wavelength region, but does not repro- 
duce the result of the kinetic equation at short wave- 
lengths even in the nonrelativistic limit. This indicates 
that the second-order terms other than the time deriva- 
tive of the dissipation terms are not useful for physical 
description of relativistic dissipative fluid. 
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Appendix A. THE ANALYTICAL CONTINUA- 
TION 

Appendix A. 1. THE ANDERSON AND WITTING 
MODEL 



In Sec. 12.31 the integrand of the dispersion relations 
Eq. (|68] | has poles. As a result, when the decay rate 
-Im a> grows to be more than l/r, we must add correc- 
tion terms of the analytical continuation. In this section, 
we derive these terms. 

For simplicity, we use the transverse shear flow mode 
(6ux - 0, 5m_l 7^ 0). The shear flow mode is given by 



(A.l) 



d P ± 

^p p 



I - ia> + /k ■ -ir 



6n 



K' 



M .6T 
-1 +z+—C\ — 



+ (y^MA + y^M± 



dp ± 
l^P P 



fo 



I — ia> + /k ■ -ir 



p^5u_i_. 



We use cylindrical coordinates in momentum space, set- 
ting k as the polar axis, and we assume that k is parallel 
to the X-axis. If we set p^ - sjip^')^ + {p~)^, the volume 
element in the momentum space (Pp reduces to 



(f'p - p±dp_i_d(pdp'\ 



(A.2) 



where is the angular coordinate around the pole k. 
We integrate with respect to and Eq. dA.U reduces 

to 

riT r°° r°° , , p'^e-^ 
—pj— dp I dp^p^- ^———xp^6uy 

4^-^X2 J_oo Jo l-ia> + ikp'-'/z 
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This integrand has first-order poles on the real axis when 
Im w < - 1 ; this singularity is removed to the upper half 
of the complex plane when Im w > -1 due to the col- 
lision term in the Anderson- Witting model Eq. $55[ . As 
in the Landau method, we consider this problem as the 
Cauchy problem. Therefore, we analytically continue 
the integrand and distort the momentum contour below 
iz(l - iciJ)/k in the complex plane. 

We have to obtain the residue of the integrand of 
Eq. ( IA.3I 1. If the integrand is fiizj/fiiz) and the pole 
is at zo, we can use the formula 



Mzo) 



f2iZ0) + 0, 



(A.6) 



where a_i is the residue. 
In Eq. (IA.3b . is given by 
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(A.7) 



Since we consider poles, the denominator /2(zo) is 0: 



ibz = 0. Then, the above equation reduces to 
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Using this result, we calculate the correction term of 
Eq. (lA3T l: 
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Using this equation, we can replace the integral variable 
p'^ with z- Then, Eq. (IA.3b reduces to 
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We integrate Eq. ( IA.12l i and multiply by 2K2k. Then, 
the equation reduces to 



-^e-'' [(1 + b^){(f + A(P + \2(f + 24-d + 24) 



- d^(d^ + 2d + 2)] 6uy, 



d^ 



(A. 13) 



(A. 14) 



This is the correction term of the analytical continuation 
ofEq. dSOll. 

In the same way, we can obtain the correction terms 
of Eqs. ( iTOt . ( l76b . and ( l82l i. First, the correction term of 
Eq. ( TTOt multiplied by kK2 is given by 
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Sscov^d, the correction term of Eq. (l76T l multiplied by 

K2 is given by 
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Finally, the correction term of Eq. (|82] | multiplied by 
kK2 is given by 
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Appendix A.2. THE MODIFIED MARLE MODEL 

In Sec. 12.21 the dispersion relations of the modified 
Marie model Eqs. ( l33l l. (l40b . ( l44b . and ( |46] | include cor- 
rection terms of analytical continuation. In this section, 
we derive these terms. 



For simplicity, we use the transverse shear flow mode 
(6x - 0, 5m J. 0). The shear flow mode is given by 



n6u± 
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We use cylindrical coordinates in momentum space, set- 
ting k as the polar axis, and we assume that k is parallel 
to the X-axis. Setting p^ - -\J(p^')^ + (p^)^, the volume 
element in the momentum space d^p reduces to 



dp — p±dp±d<pdp'^ 



(A.19) 



where (p is the angular coordinate around the pole k. 
We integrate with respect to 0, and Eq. ( lA.lb reduces 

to 
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As in the Anderson- Witting model, we consider this 
problem as the Cauchy problem. Therefore, we analyt- 
ically continue the integrand and distort the momentum 
contour below {u)z + 1^X2! Ki) I k. 

We have to obtain the residue of the integrand of 
Eq. ( IA.2OI 1. As in Sec Appendix A. 1 we use Eq. ( IA.6I 1. 
In Eq. ( IA.20b . /j is given by 
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Since we consider the pole, the denominator /2(zo) is 0: 



(A.24) 
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From the scalar product of the four-momentum p^, we 
obtain 

(A.25) 



Using the above equations, we get 




(A.26) 



(A.27) 



where we determine the sign of p''^ and z to be z ^ oo 
and Im p'* < when p'j^ — > oo. 

Using these equations, we can calculate correction 
terms of analytical continuation. Considering that 
we have to analytically continue when p''-' becomes 
Im p'" < 0, the correction terms of Eq. ( IA.20l i are given 
by 



no 
HK, Jo 

X 6d {p 

no 2n 



Jo J-c< 



'2.711 

dp'-'—6(-lm p'-') (A.28) 
ikz 



k z 



4^K 



P'l e 



2.71 

dp'Ji-lmp'- 

1 K Jo 



z 1 



01 £^ 

k z 



6uv 



where 6(x) is the Heaviside step function given by 



eix) =0 if x<Q, 
e{x) =1 if x>Q. 



(A.29) 
(A.30) 



Using Eq. ( IA.28b . we replace the integral variables p'^ 
with z- Then, the correction terms are 
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4^Ki T 
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(A.31) 

where c is the value of z evaluated when Im p'^ip'^ - 
0. We integrate of Eq. ( IA.31b and multiply by 2kK\ 
to obtain the correction terms of Eq. ( l44l i. In the same 
way, we can obtain correction terms of Eqs. (l33l l. (l40l l. 
and (|46] |. 



Appendix B. RELATIONSHIP BETWEEN THE 
BGK MODEL AND THE MATCH- 
ING CONDITIONS 

In general, the local equilibrium distribution func- 
tion lacks a physical meaning until it fulfills the match- 
ing conditions. Although the matching condition of the 
nonrelativistic BGK model is unique, that of the rela- 
tivistic BGK model has different forms depending on 
the fluid four-velocity. For example, the matching con- 
dition of the Marie model Eq. ( fTTI l diff'ers from that of 
the Anderson-Witting model Eq. ( l56l l. In this section, 
we consider this diff'erence. 

The kinetic equation has to fulfill the conservation 
laws of particle four-flow, the energy-momentum ten- 
sor. When we consider the kinetic equation of the BGK 
model, the conservation laws reduce to 

^"^"^ / ^^^"^^ / 7re(/>A) = o, (B.i) 



j p' 
J p' 
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'^^P'Q(fJe,)-0, 



where Q{f, feq) is the collision term of the BGK model. 
The collision term Q depends on the local distribution 
function, so the above equations become constraints on 
it. Moreover, the number of conservation laws is the 
same as the number of degrees of freedom of the lo- 
cal equilibrium distribution function. As a result, the 
matching condition of the BGK model is determined by 
the conservation laws. 

First, we consider the matching condition of the 
Marie model. In this case, the above equations reduce 
to 



P° Tm 



1 



0, (B.3) 



- r ^P"^ = -- (M' - O = 0, (B.4) 



where n is the number density of particles, and e is the 
energy density per particle. Eqs. ( IB.3l l and (IB.4l i indi- 
cate that in the Marie model, the conservation law of the 
particle four-flow Eq. ( l33l l becomes the matching con- 
dition of n/e, and the conservation laws of the energy- 
momentum tensor Eqs. (l40b . (l44l i. and (l46l l become the 
matching condition of particle number density n and the 
Eckart velocity u . Note that Eq. ( IB. 3b is not strictly 
the matching condition of the energy, so that in general 
e - ^ Oin the Marie model Q\. 
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In the nonrelativistic case, the conservation laws of 
particle four-flow are equivalent to the matching condi- 
tion of the particle density n, and that of momentum flux 
is equivalent to the matching condition of the fluid ve- 
locity u. Thus, from these two equations we can derive 
a continuous equation: -ia)6n/n + /k ■ 6u = 0. In the 
Marie model, however, we have to use all the equations 
to derive a continuous equation; Eq. (l33T l xi(/Kik+ 
Eq. (l40l i - Eq. (l46b xco^jk reduces to 



6n 

Kii^oj ^Ki^-5\x 

n 



0. 



(B.5) 



The above equation is equivalent to the continuous 
equation. 

Next, we consider the Anderson- Witting model. As 
in the Marie model, we substitute Q for the collision 
term in the Anderson-Witting model. Then, the conser- 
vation laws reduce to 



( 



(B.6) 



(B.7) 

Eqs. ( IB. 6b and ( IB. 71 ) indicate that in the Anderson- 
Witting model, the conservation law of the particle four- 
flow Eq. (iTOt becomes the matching condition of parti- 
cle density n, and the conservation law of the energy- 
momentum tensor Eqs. ( l76l l. dSOb . and ( l82b becomes 
the matching condition of energy density e and the 
Landau-Lifshitz velocity u . In the Anderson- Witting 
model, the velocity is not the Eckart velocity but the 
Landau-Lifshitz velocity. Therefore, we cannot derive 
the continuous equation from Eqs. ( fTOl i. ( l76l l. (iSOl l. and 
( l82b but we can derive the conservation law of energy 
-i(ji)6(ne) + inhk ■ 6u - 0, where h is the enthalpy per 
particle of the unperturbed gas. To obtain this equation, 
we use only the conservation laws of energy and mo- 
mentum; Eq. ( |76] | X/-I- Eq. ( |82] | xb reduces to 



6n 

- iaj(3K2 + (Ki) — 
n 



(B.8) 



+ i^K^k ■ 5u = 0. 



6T 



After some calculation, we can check that the above 
equation is equivalent to the energy conservation laws: 



Appendix C. ASYMPTOTIC ANALYSIS IN THE 
LONG WAVE LENGTH LIMIT 

To study the behavior of the roots of the dispersion 
relation, we take the long wavelength limit. For sim- 
plicity, we analyze the Anderson- Witting model in this 
section. First, we consider the shear flow mode. From 
Sec. 12.31 the dispersion relation of the shear flow is 



J A 

J P" 



(C.l) 



We expand the integrand in powers of -io) + ikv^ and 
neglect terms higher than second order on the right-hand 
side. Then, the above equation reduces to 



n6ux_ = J -^Pi-\ 



1 + {ioj - ikvx) 



(C.2) 



+ {-a? - l^vl + 2ukv,)]f,, 



X 



6n I ^ nO K' \ST 



Rewriting the above equation yields 

73r 



= 1^(1 + ico — a/) — ak?^ n6u^, 
C d^P 2 If 

Neglecting a?, we obtain 

CD — -iak^. 



(C.3) 



(C.4) 



(C.5) 



This reproduces the dispersion relation of the shear flow 
mode in the long wavelength limit. 

Next, we study the long wavelength limit of the ther- 
mal conduction and sound wave modes. From Sec. 12.31 
the conservation of particle number Eg.dTOli is 



6n r d^p Q 



feq 



ia)5{ne) + inhk. ■ 5\x- Q. 



(B.9) 



ia> + ikVf 
6n I , K'^ \ST 



(C.6) 
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As in the shear flow mode, we expand flie integrand in 
powers of -iu) + ikv^ and neglect terms higher than sec- 
ond order on the right-hand side. Then, the above equa- 
tion reduces to 



n J 

+ {-co" - k\l + 2ukv,)]feg 



(C.7) 



6n I , / K'^ \6T 



5n 



= (1 -I- iu - - {vi)k^)— + i-ikT" + 2cokT'")6u, 
n 



6T 



eg- 



(C.8) 



where 



A^c,nTr^\ (C.13) 

B = [2c„«rr"(-r°" + + sr") (c. 14) 

+ t'^\t"" + (-3 + c,nT)T'' - Ac,nT{vl)) 



K2 



{-2c,nT(T-"f + T°\T" + c,nT{x^,)))z 



K2 



(C.15) 



r0xx/^,2\ 



(C.16) 



We obtain three roots, denoted by a)T,(Ds± to second 
order in k: 



Rewriting the above equation yields 



6n 



(ioj -Lo^- <v^)r)— -I- i-i + 2w)/tr"(5M^ 
n 



(C.9) 



+ k' 



ST 



Similarly, the conservation of energy Eq. ( l82l i reduces 



to 



[(io) - a?)c„nT - T"k^] — + (-/ + 2oj)knh6u^ 
n 



(C.IO) 



+ [(iw - a/)nTcy 



+ k'T"\'i + YC\-k'T 



T 



where n, T, and h are the particle number density, tem- 
perature, and enthalpy of the unperturbed state, respec- 
tively, and Cy is the heat capacity per particle. 

From Sec. 12.31 the dispersion relation of the 
Anderson- Witting model includes the conservation of 
energy Eq. ( IB.9I ), so we use the conservation of energy 
instead of the conservation of momentum Eq. ( |76] |. 

- iio{CyT6n + nCv6T) + inhkSux — 0- (C.l 1) 

From the above equations, we can obtain the dispersion 
relation in the form 



iAkcs" + Bk^(J + iCk^u - Dk^ = 0, (C. 12) 



D , 

ut = -i-^^, (C.17) 
/C AD-BC , 

The first solution represents the thermal conduction 
mode, and the second two represent the sound wave 
mode. 
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